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The Fermi-hAT view of cosmic rays and interstellar gas in the Cygnus region: a not so special 
spot of the Local Arm 
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Abstract: The Cygnus region hosts the most conspicuous star-forming region close to the Sun, embedded in a giant 
complex of molecular clouds in the Local Arm. We present an analysis of the Fermi-hAT observations of Cygnus 
intended to probe the cosmic-ray and interstellar-matter content of the region. From gamma-ray data we estimate a total 
of 8tf X 10^ Mq of interstellar gas in the complex at a distance of 1.4 kpc. The gamma-ray emission from the atomic gas 
supports the average H i spin temperature derived from radio absorption/emission pairs to estimate its column densities. 
The Xco ~ N{}i2)/Wco ratio derived in the massive Cygnus complex is consistent with other LAT estimates for 
clouds in the Local and Perseus arms. The mass of dark gas, escaping H i and CO observations but traced by dust and 
gamma rays, amounts to ~ 10% of the total. We find an average gamma-ray emissivity per interstellar H atom in the 
0.1-100 GeV energy band in good agreement with measurements in other segments of the Local Arm. We infer that the 
cosmic-ray population averaged over a few hundred parsecs is fairly uniform in density and spectrum along the Local 
Arm. Despite the presence of potential accelerators and much larger interstellar densities in Cygnus compared to the solar 
neighborhood, their cosmic-ray populations are similar on such a scale. 
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1 Introduction 

Interstellar 7-ray emission produced by cosmic-ray (CR) 
interactions with interstellar gas, through nucleon-nucleon 
collisions and electron Bremmstrahlung, is a tracer of the 
product of CR and interstellar-matter densities throughout 
the Galaxy. The data which are being collected by the 
Fermi Large Area Telescope (LAT) |6| are therefore pro- 
viding a wealth of new information on the interstellar envi- 
ronment of the Milky Way. 

We used LAT data to characterize the CR and gas content 
of the Cygnus region, which is located around the Galactic 
longitude I — 80°, tangent to the Local Arm. A super- 
massive molecular complex harbors the region of massive 
star formation known as Cygnus X, at ^ 1.4 kpc from the 
solar system ifTTl . Beyond the Local Arm two segments of 
the Perseus and outer spiral arms are seen in this direction. 

In this contribution we present the global model that was 
sought for the high-energy 7-ray emission from Cygnus 
and we discuss the implications for the properties of CRs 
and the census of the interstellar medium (ISM) over the 
scale of the whole complex, i.e. a few hundred parsecs. 
The innermost Cygnus X region will be discussed in an- 
other contribution to this conference 1,22]. 



2 LAT Data Analysis 
2.1 7-ray Data 

We used the first two years of the LAT survey applying 
stringent event selection criteria to limit the background 
due to misclassified CR interactions in the LAT |4| and 
the contamination from the Earth's atmospheric emission. 
We analysed LAT data in the region at Galactic longitude 
72° <l < 88° and latitude |6| < 15°. Photons were used 
to construct maps on a 1/8° angular grid and over ten bins 
in energy from 100 MeV to 100 GeV. 

The 7-ray emission from this region is dominated at en- 
ergies < 10 GeV by three bright pulsars: J2021H-3651, 
J2021H-4026 and J2032H-4127. To increase the sensitivity 
to diffuse and fainter extended sources at these energies for 
each pulsar we removed from our data sample the photons 
detected in the periodic time intervals where the pulsed 
emission peaks, in a region around the source position 
with a radius approximately equivalent to the acceptance- 
averaged 95% point-spread function (PSF) of the LAT as a 
function of energy. 
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2.2 Analysis Model 

Since the Galaxy is transparent to high-energy 7-rays and 
the CR densities are thought to be uniform over the charac- 
teristic dimensions of interstellar complexes, we can model 
to first order emission from gas as a linear combination of 
gas column densities for the different complexes along the 
line of sight and the different phases of the ISM. 

Column densities of atomic hydrogen, Hi, were derived 
from the observations of its 21 -cm line lfT2l[T3]| . Column 
densities of molecular gas are assumed to be proportional 
to the integrated intensities of the 2.6-mm line, Wco |7 8 1, 
(H2 ) — Xqo ■ Wco ■ H I and CO data were used to build 
maps of gas column densities partitioned along the line of 
sight to separate the Cygnus complex from the spiral arms 
in the outer Galaxy, according to the method described in 

isia. 

The dark neutral gas which is missed by H i and CO obser- 
vations was recovered by a dust excess map {\V\, namely 
by the visual extinction Ay estimated from the 2MASS 
source Catalog |[10l[19l|20l minus the best-fit Unear combi- 
nation of the afore described H i and CO maps. 

The emission from interstellar gas is therefore modeled as 
a linear combination of the Hi, CO and Ay.cxc maps with 
free combination coefficients independently in each energy 
bin. These combination coefficients are the emissivities, 
i.e. 7-ray emission rates, per H atom, qui, per Wco unit, 
gco, and per Ay ,cxc unit, (7Av ■ 

The interstellar emission model is completed by the large- 
scale Galactic inverse-Compton emission produced by 
CR electrons interacting with low-energy radiation fields, 
which was modeled using the GALPROP codfl run 
54_87Xexph7S. An isotropic background component, 
with free normalization in each independent energy bin, 
was added to take into account the isotropic 7-ray emission 
and the residual background from misclassified CR inter- 
actions in the LAT. 

Several sources are included in the model: all the sources in 
the first LAT source Catalog (IFGL) which are identified, 
associated at high confidence level with Active Galactic 
Nuclei or show strong variability. We then introduced some 
extended sources detected with LAT data: the supernova 
remnants known as the Cygnus Loop lfT4l and 7 Cygni, 
the source VER J2019H-407 fTS\ and the extended source 
in the innermost part of the Cygnus X region discussed 
in II22II . The remainder of the three bright pulsars, phase- 
suppressed as described in § 12.11 was modeled for each of 
them by a point source to take into account the off-pulse 
emission and by a second source, for which, after folding 
with the instrument response as described in §|23] expected 
counts were set to null within the exclusion radius intro- 
duced in ij |2.1l to take into account the spill-over in the PSF 
tails. For each source the flux was let as a free parameter 
independently in each energy bin. 



2.3 Analysis Method 

The standard LAT analysis softwar^ was used to calculate 
the energy-dependent exposure, which was then scaled to 
take into account the selection on pulsar phases (i; 12. 11 1. The 
analysis model described above was corrected for the expo- 
sure and convolved with the PSF, to obtain the number of 
counts expected from the model in each pixel of the maps. 
This was then fit to the LAT data by a maximum-likelihood 
method based on Poisson statistics in order to determine 
the best model parameters. 

The analysis was performed using the Instrument Response 
Functions (IRFs) of the series P6_V3 [18|. The relative 
systematic uncertainties described therein are taken into ac- 
count in the interpretation of the results. 

3 Results and Discussion 

3.1 The Issue of the H i Spin Temperature 

The preparation of the A^(Hi) maps described in i) l2.2l re- 
quires to account for the optical depth of the gas, which is 
usually done by assuming a uniform value of the so-called 
Hi spin temperature, Ts- No simple solutions of the Hi- 
line transfert equation exist beyond the approximation of 
uniform Ts- We therefore explored a wide range of possi- 
ble average temperatures from 100 K to 00 (optically thin 
medium). The corresponding maximum-likelihood values 
obtained from the 7-ray analysis are shown in Fig. [T] 
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Figure 1: PRELIMINARY - Maximum likelihood ob- 
tained by fitting the model to the LAT data for some dis- 
crete values of the H i spin temperature Ts used to prepare 
A^(H i) maps: 100 K, 125 K, 250 K, 400 K and the Hmit of 
infinitely high T5 (optically thin medium). 

The results support the average Ts values of a few hun- 
dred K estimated from radio absorption-emission pairs 19] 
to estimate the H i column densities. The spin temperature, 

1. http :/ /galprop . Stanford . edu/ 

2. http : // f ermi .gsfc.nasa. gov/ ssc/data/ 
analysis/ documentation/ Cicerone/ 



32nd International Cosmic Ray Conference, Beijing 201 1 



however, is known to be not uniform: we used the differ- 
ent values explored to gauge the subsequent uncertainties 
affecting the results. 

3.2 CR Densities in Cygnus and along the Local 
Arm 

Fig. |2] shows the spectrum of 7-ray emissivity per H atom 
in the Cygnus complex, which is consistent with the pre- 
dictions based on the CR spectra directly measured near 
the Earth and 7-ray emission from the nearby interstellar 
space as observed by the LAT fT\. 
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Figure 2: PRELIMINARY - Spectrum of emissivity per H 
atom in Cygnus. Shaded bands indicate the systematic un- 
certainties, combining the imperfect knowledge of the in- 
strument response and of the H i optical depth. The dashed 
line shows the expectations based on a CR spectrum con- 
sistent with direct measurements near the Earth [IJ. 

Fig. |3] compares the Hi emissivity integrated above 
200 MeV in the Cygnus complex and other segments of the 
Local Arm and regions toward the outer Galaxy lf3ll5l fT5l . 
Given the large uncertainties due to the issue of the H i spin 
temperature (§[TTJ, Fig.|3]shows that the CR population av- 
eraged over a scale of a few hundred pc is uniform within 
10% - 35% along the Local Arm. 

3.3 The Census of the Interstellar Medium 

By the definition A^(H2) = Xqo • T^co. we obtain 
qB_2 = ^co ' Ico- Assuming that the CR densities are 
the same in all the phases of an interstellar complex, we 
expect = 2gHi- We therefore expect the emissivity per 
Wco unit to linearly scale as the emissivity per H atom, 
qco — 2Xco • Qui- This is verified in Fig. |4] (left), with 
good confidence in the 0.1-10 GeV energy range. A small 
(< 3a) deviation is seen > 10 GeV (i.e. at low emissivi- 
ties), possibly due to the hard extended emission Il22l ob- 
served in the innermost Cygnus X region where most of the 
CO-bright gas is located. From the slope of the best-fit lin- 
ear relation we can derive Xqo and estimate the masses of 
CO-bright gas. 
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Figure 3: PRELIMINARY - Emissivity per H atom inte- 
grated for E > 200 MeV as a function of Galactocentric 
radius, combining results from this analysis (red, or light 
gray) and from the studies of the second |I3] black] and third 
1 5 , blue or dark gray] Galactic quadrants. Statistical errors 
are comparable with marker dimensions, whereas the dif- 
ferent markers show the results for different values of the 
H I spin temperature Ts used to prepare (H i) maps. 

We similarly expect the emissivity per Ay.oxc unit to lin- 
early scale as the emissivity per H atom, as we verified 
in Fig. |4] (right). The linear correlation implies that 7-ray 
emission associated to dust excesses and atomic gas is pro- 
duced by the same physical processes, i.e. that Ay.oxc in- 
deed trace interstellar gas. From the best-fit linear rela- 
tion we can derive the dust-to-gas ratio in the dark neu- 
tral phase, = ^(H)/Av,exc, and estimate its mass. 
Whether the dark neutral gas is CO-quiet molecular gas 
or atomic gas missed because of the approximations in the 
handling of the radiative transfert of the H i line should be 
investigated in the future. 

Assuming a distance of 1.4 kpc from the solar system, the 
total masses estimated on the basis on the LAT data anal- 
ysis amount to 8lf x IO^Mq in the Cygnus complex, of 
which ^ 60% in the form of atomic gas, ^ 30% in the form 
of CO-bright molecular gas and ^ 10% in the form of dark 
neutral gas. 

The Xqo ratio obtained for the Cygnus complex appears 
consistent with estimates from LAT data in other segments 
of the Local Arm and in the Galactic disc over a few kpc 
beyond the solar circle (Fig|5]l. It is significantly larger than 
the measurement for the nearby clouds of Cassiopeia and 
Cepheus |l3]; the origin of this difference will be investi- 
gated in the future. 

4 Summary 

The Cygnus region appears as a very peculiar segment of 
the Local Arm because of the large masses of interstel- 
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Figure 4: PRELIMINARY - Emissivity per Woo unit (left) and per Ay.cxc unit (right) versus emissivity per H atom for 
the different energy bins considered in the analysis. The dashed lines represent the best-fit linear relations taking into 
account uncertainties on both axes. 
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Figure 5 : PRELIMINARY - Xco as a function of Galacto- 
centric radius. Points represent LAT measurements for the 
Cygnus clouds (red or light gray) and several complexes 
seen in the second [ 3 , black] and third |T blue or dark gray] 
Galactic quadrants. Some models reported in the literature 
llT6ll2n are shown for comparison. 



lar gas, confirmed by its 7-ray emission, and the intense 
massive-star formation activity, and therefore the potential 
presence of numerous particle accelerators. Despite these 
peculiarities, the properties of CRs over the scale of the 
whole complex are similar to those of other segments of 
the Local Arm studied with LAT data. 
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